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Abstract 

The experimentally observed temperature dependence of the quartic coeffi- 
cients in the Landau-Devonshire expansion for BaTi03 is naturally accounted 
for within a proper fluctuation model. It is explained, in particular, why one 
of the quartic coefficients varies with temperature above T c while the second 
is constant. The tetragonal phase in BaTi03 is argued to exist essentially due 
to the thermal fluctuations, while the true Landau-Devonshire expansion with 
temperature-independent coefficients favours the rhombohedral ferroelectric 
phase. Certain conclusions concerning the temperature dependence of the 
sextic Landau-Devonshire coefficients are also made. 

PACS numbers: 77.80. Bh, 77.84.Bw 
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The phase diagram of the barium titanate contains three lines of phase transitions and its 
structure is known to be properly reproduced by the phenomenological Landau-Devonshire 
theory jl] . It is accepted within the phenomenological approach that all the coefficients in the 
Landau free energy expansion should be either constant or weak functions of temperature, 
pressure, etc., apart from the quadratic term that changes sign crossing the second-order 
transition line or the low-temperature spinodal. For the barium titanate, however, the 
matching with the theory of the experimentally observed temperature dependences of the 
nonlinear dielectric susceptibility and spontaneous polarization forces one to allow for the 
strong temperature dependence of some higher-order coefficients Being in conflict 

with the spirit of the Landau theory itself, rapid temperature variation of the quartic and 
sextic coefficients is also quite unexpected from the microscopic point of view. Indeed, the 
barium titanate is a displacive ferroelectric with weak anharmonicity that can only results in 
rather slow temperature dependence of the macroscopic parameters ||. Moreover, there is an 
extra point showing how unsatisfactory the real situation is: because of the fast temperature 
variation of the certain coefficient, the six-order Landau-Devonshire expansion for BaTiC>3 
with updraded coefficients |||7|] turns out to loose its global stability at T = 443 K, i. e. at 
the temperature that exceeds the cubic-tetragonal transition point T c by only 50 K. 

In this Letter, we will show that the temperature dependence of the quartic Landau- 
Devonshire coefficients, experimentally observed in a paraelectric phase of BaTiC>3, can be 
attributed to the thermal fluctuations of polarization and naturally accounted for within 
a proper fluctuation model. It will be explained, in particular, why in experiments one 
of the quartic coefficients demonstrates well-pronounced variation with temperature while 
the second is temperature independent. The conjecture will be put forward concerning the 
structure of the true Landau-Devonshire form for BaTiC>3 with the constant (temperature 
independent) coefficients, and the first-order fluctuation corrections to the sextic Landau- 
Devonshire coefficients for T > T c will be found. 

Barium titanate is a displacive ferroelectric undergoing the first-order ferroelectric phase 
transitions, with the six-order anharmonicity and electrostriction playing an essential role in 
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forming its phase diagram. The electrostriction is known to convert the second-order transi- 
tion, appropriate to the clamped crystal, into the first-order one. Strong dipole-dipole inter- 
action affects the vibrational spectrum of BaTiOs, resulting in a big gap between transversal 
(soft) and longitudinal polarization modes. These main features are properly described by 
the effective Hamiltonian that is the natural generalization of the Landau-Devonshire free- 
energy expansion: 
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H int = J dPx^Mt + 4 + 4) + lMl<t>l + 0203 + 030?) + 7l(01 + 02 + 03) 

+ 3 7 2 (0^ + 44 + 0203 + 0203 + 0301 + 44) + ^3444} , (2) 

H str = J rf 3 x{gn(M n 0? + W 22 02 + «3303) + 5l2[«ll (02 + 4) + U 22(4 + 4) 

+ U 33 (4 + 4)} + <?44(«230203 + %10103 + «12020l)}- (3) 

Here <p a and u a p are the Cartesian components of the fluctuating polarization and strain, <p aq 
stands for the Fourier transform of (f) a (x), n a = ^, C is the Curie constant, A 2 ~ e 1 ~ sa^ 1 
is the dipole gap in the fluctuation spectrum, a being a lattice constant. Contrary to the 
original Landau-Devonshire expansion, the Hamiltonian (1-3) accounts for the inhomoge- 
neous fluctuations of the polarization and elastic strains making it possible to explore the 
fluctuation effects in BaTi03. 

Dealing with the first-order phase transition, we are in a position to study the fluctuation 
effects in the region where thermal fluctuations of the order parameter are weak. Hence, in 
what follows we limit ourselves by the calculation of the first-order fluctuation corrections to 
the quantities of interest. The quantities to be found are the full four-leg and six-leg vertices, 
reducing, under zero external momenta, to the effective ("dressed") quartic Bi and sextic 
Ti coefficients, entering the Landau-Devonshire expansion. Among five Feynman diagrams 
representing the lowest-order non-trivial terms in the perturbative expansions of the vertices 



mentioned, only two one- loop graphs give contributions that rapidly, as (T — To) -1 / 2 , grow 
up approaching the transition point. These graphs have an obvious structure and may be 
calculated in a standart way using the propagator 

r / \ _ k B T(5 aP - n a np) 

GaM ~ (Ce )-i(T-T ) + s^ (4) 

with a longitudinal part fully neglected. The results for the quartic couplings B\ and B 2 are 
found to be: 



B 2 = fa- ^J^T &fi + 18 ^ 2 + 13 ^' ( 5 ) 

The polynomials in brackets are easily seen to coincide with those of the one-loop contri- 
butions to the renormalization-group /5-functions of the cubic ferroelectric ||||. It is not 
surprising since, in fact, the same integrals and tensor convolutions are evaluated in both 
cases. The fluctuation correction to B\ given by the first equation (5) is consistent with the 
results of Vaks @,fL(J ? who first evaluated this correction and showed that it is essential for 



the case of BaTiC^. 

To proceed further, we have to estimate the "bare" coupling constants (3\ and P2 for 
BaTi03 trusting upon the experimental data available. Aiming to extract the necessary 
information from experiments, one should realize that i) what is measured are not bare but 
dressed couplings with the fluctuation contributions included, and ii) in experiments, the 
Landau-Devonshire coefficients are measured for stress-free (not clamped) crystals. Hence, 
in order to estimate pi and fa, we have first to express them via their analogs for free crystal, 
Pi and P2 ■ This problem is solved by evaluation of the elastic strains caused by the non-zero 
uniform polarization and consequent renormalization of coefficients in the relevant Landau- 
Devonshire expansion Jl|. All the elastic and electrostrictive moduli are known for BaTiC>3 
0, making corresponding calculations straightforward. The final result is as follows: 

Pi = p{ + 7.4 • 10 8 Vm 5 C- 3 , p 2 = Pi - 2.3 • lOVm 5 ^ 3 , (6) 



Since in the vicinity of T c the elastic and electrostrictive moduli weakly depend on tem- 
perature, similar relations should be valid for the fluctuation modified (dressed) quartic 
coefficients B\, B 2 , B(, and B 2 . 

At the transition point (T c = 393 K), B[ = -2.0 • 10 8 Vm 5 C~ 3 , B{ = 1.6 • 10 8 Vm 5 C- 3 
§J7|] and, hence, B l = 5.4 • 10 8 Vm 5 C- 3 , B 2 = -0.7 • WVnt'C' 3 . It is easy to see that for 
B\ > B2 the rhombohedral phase has lower free energy than the tetragonal one, provided 
the six-order form is isotropic, i. e. does not influence their competition. Hence, in the 
clamped crystal the quartic form of the Landau-Devonshire expansion strongly favours the 
phase transition into the rhombohedral phase. The same is true for higher temperatures 
T = 415 K and T = 423 K, where B{ = -1.3 • 10 8 Vm 5 C- 3 0, B l = 6.1 • 10 8 Vm 5 C- 3 
and B( = -1.0 • 10 8 Vm 5 C- 3 [§, B 1 = 6.4 • 10 8 Vm 5 C~ 3 , respectively, with B{ and B 2 kept 
unchanged. 

Now we are ready to estimate to what extend the fluctuations can modify the behaviour of 
barium titanate in the vicinity of T c . As is seen from Eqs.([5]), the quartic Landau-Devonshire 
coefficients should vary with temperature in a similar way, provided the fluctuation correc- 
tions to them are of the same order of magnitude. Let us compare the magnitudes of the 
fluctuation terms B^p and B^ within the domain where parameters (3\ and j3 2 have relevant 
values and signs. The quantity characterizing relative weights of B% and B^ is their ratio, 

that can be found directly from Eqs. (|5|). The plot of the ratio R> > = — as a function of 

a 1 b (1) 

r = j± is shown in Fig. 1. Analyzing this function, one can find, in particular, that -^pj 

does not exceed 0.1 if I 1 lies between —28 and —2.3. 

P2 

Let us estimate further |^ in BaTiOa. Since, according to experiments, B\ varies with 
the temperature appreciably, an estimate for (3\ is expected to have somewhat limited ac- 
curacy. To keep the perturbation theory more or less meaningful, one should adopt that 
the fluctuation term B^ is, at least, two times smaller than the value of B\ at T = T c . 
At the transition temperature B\ = 5.4 • 10 8 Vm 5 C~ 3 and this coefficient was shown to 
increase when the temperature grows up. Hence, the limitation Pi < 10 9 Vm 5 C~ 3 looks 
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quite reasonable. Another quartic coefficient, B 2 , does not depend on T and the estimate 
/?2 — B 2 = —0.7 ■ 10 8 Vm 5 C~ 3 may be considered as an accurate one. Since, in any case, 
Px > B^T = 393K) = 5.4 ■ 10 8 Vm 5 C- 3 , we see that -14 < j| < -8 in barium titanate. 
It implies that, according to Fig. 1, B 2 ^ < 0.083^, and the fluctuation correction B 2 ^ is 
inevitably very small. This means that B 2 practically does not depend on T and explanes 
why in BaTiC>3 the Landau-Devonshire coefficient B 2 is insensitive to the temperature. 

Apart from the smallness of the theory naturally accounts for the experimentally 
observed sign of the fluctuation contribution to B\. Indeed, as seen from Eqs. (fjD, a sign of 
b[^ is completely controlled by the first term proportional to (3f: big number 24 provides 
a positiveness of the polynomial in brackets avoiding any possibility for the rest terms to 
compete with the first one under any values of f3\ and f3 2 . Hence, the theory definitely 
predicts that the fluctuation correction b[^ is negative and, therefore, approaching T c from 
above B\ and B( should vary downward. This conclusion is in agreement with experiments 

As we have already seen, the experimental data available do not allow to fix the true, 
temperature independent value of the Landau-Devonshire coefficient f3\. However, there 
are obvious requirements that enable us to improve the crude estimate for f3\ presented 
above. The theory developed, based on the first-order perturbative calculations, is believed 
to account for the main features of the behaviour of BaTiC^. Hence, it should explain 
the variation of B{(T) by, at least, 10 8 Vm 5 C~ 3 [[J and remain valid, at worst, at the 
semiquantitative level. It is possible to meet both requirements only accepting that the 
true value of (3\ is appreciably bigger than 6.4 • l0 8 Vm 5 C~ 3 and appreciably smaller than 
1.1 • W 9 Vm 5 C~ 3 , i. e. it lies somewhere between 7 • 10 8 Vm 5 C~ 3 and 10 9 Vm 5 C~ 3 . As a 
result, the stress-free crystal should possess the Landau-Devonshire coefficient that obeys 
the inequalities -0.4 ■ 10 8 Vm 5 C- 3 < p[ < 2.6 • 10 8 Vm 5 C- 3 . 

The upper part of this interval is of particular interest. The point is that, whenever (3( 
exceeds 1.6 • 10 8 Vm 5 C~ 3 = /3 2 , the crystal with the fluctuations being "switched off" would 
undergo a phase transition into the rhombohedral phase rather than into the tetragonal one. 
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If this were true, i. e. the inequality (3{ > {5 2 took place, the tetragonal and orthorhombic 
phases in barium titanate would exist essentially due to the thermal fluctuations, while the 
rhombohedral phase would survive only at sufficiently low temperatures where the fluctua- 
tions are weak enough. Thus, the analysis based upon the effective Hamiltonian (1) turns 
out to support the conjecture about the fluctuation stabilization of the high-temperature 
ferroelectric phases in BaTi03 formulated first within the microscopic theory ||. 

This conjecture looks rather attractive. It is worthy to discuss it in more detail. It turns 
out that, apart from the results the of first-order calculations, there exist two extra arguments 
in favour of the scenario just described. The first is as follows. As we have already seen, 
the lowest-order correction to j3{ is considerable and, therefore, the higher-order fluctuation 
contributions can influence the results appreciably. The leading perturbative term shifting 
the first-order estimates is the second-order one and it has a sign opposite to that of B[ . 
One can show that account for the positive second-order fluctuation term in the course of fit 
of the experimental data results in the shift of the estimated value of (3( upward. Hence, 
within the refined theory the fluctuation origin of the tetragonal and orthorhombic phases 
in BaTiC>3 would become more plausible. 

The second argument deals with the fluctuation corrections to the six-order coefficients 
in the Landau-Devonshire expansion. Let us calculate them. Evaluating the six-leg one-loop 
diagram, we arrive to the following expressions for r,: 

Ti = 7i - /(T)(120/? l7l + 6/3 l72 + IO&71 + I8&72), 

T 2 = 72 - /(T)(5/3 lTl + 63/3 l72 + 2/5 l73 + 15&7i + 63&72 + 6/^73), 

T 3 = 7s - /(T) (18/3i7a + 24/3 l73 + 54& 72 + 58/3 273 ), (7) 

where f(T) = k B T^Ce^{l0ns 3 VT - T )~\ 

Above formulas are written down for the physical case when the polarization vector has 
three Cartesian components. In fact, the fluctuation terms were calculated for the generic 
model with the n-vector order parameter. It was done in order to reserve an opportunity 
for independent check of the results obtained. Indeed, as is well known, the model with 
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cubic anisotropy possesses the special symmetry property under n = 2: if one turns the field 
(0i)02) by 45° in its two-dimensional space, the coupling constants are transformed, but 
the structure of the Hamiltonian remains unchanged. This implies some exact symmetry 
relations between coupling constants fn| , that interrelate dressed Tj) and bare 7*) 
couplings and can be used to approve (or disregard) the results of perturbative calculations 
in the arbitrary order. Having made such a check up, we found that the n-vector analog of 
Eqs.(0) do obey these symmetry relations when n = 2. 

To estimate the magnitudes of fluctuation corrections rf 1 , the experimental information 
about 7j is necessary. The experiments, however, were carried out in the ordered (tetragonal) 
phase and yielded, in particular, strong dependence of T 1 on temperature. This tempera- 
ture dependence was already mentioned to be so dramatic that, being extrapolated to the 
paraelectric region, makes the Landau-Devonshire form unstable for T > T c + 50K. It is 
hardly believed therefore that the experimental data available can be used to extract more 
or less reliable estimates for 7^ 

In such a situation it is natural to analyze the general structure of the correction terms 
Vf^ in Eqs. ([?]), aiming to find some conclusions that are insensitive to the concrete values 
of 7j. Let us proceed, accepting that 71 ~ 72 ~ 73- As we have already found, in barium 
titanate l/^l << It means that the magnitudes of r| are determined, in fact, by 

the terms in Eqs. ([F]) containing f3\. The numerical coefficients before these terms are 
seen to be markedly different. The biggest one (120) stands in the expression for 
making Ti stronger dependent on temperature than two other coefficients. Moreover, the 
structure of 1^ fixes its sign. Since (3\ > 0, 1^ is negative and the fluctuations diminish 
the Landau-Devonshire coefficient T\ provided 71 > 0. The positiveness of 71 is, in its turn, 
inevitable, because this coefficient is responsible for the global stability of the system outside 
the fluctuation region. 

So, we see that approaching T c the coefficient Ti decreases more rapidly than T 2 and T 3 . 
On the other hand, as one can see, the smaller r 1; the more stable is the tetragonal phase. 
Hence, the fluctuations modify the six-order form in the Landau-Devonshire expansion in a 
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way that favours the transition into the tetragonal phase. It confirms the conjecture about 
the fluctuation stabilization of the phases lying between the cubic and rhombohedral ones 
at the phase diagram of barium titanate. 

To conclude, we have shown that the temperature dependence of the quartic Landau- 
Devonshire coefficients in a paraelectric phase of BaTi0 3 can be explained as the fluctuation 
effect. The theory developed naturally accounts for a signs of the temperature variations 
of B\ and makes it clear why the second coefficient, B 2 , is temperature independent. The 
conjecture was formulated that the true Landau-Devonshire form for BaTiC>3 with the tem- 
perature independent coefficients favours the transition from cubic into the rhombohedral 
phase, and it is the thermal fluctuations what stabilizes the tetragonal and orthorhombic 
phases and provides the space for them at the phase diagram. 
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CAPTION 

Fig. 1. The ratio of the first-order fluctuation corrections R ( > = -^jj as a function of 

§x 



the ratio of bare coupling constants r = I 1 . 
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